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Abstract. Using a detailed radiative transfer analysis, com- 
bined with an energy balance equation for the gas, we have 
performed extensive modelling of circumstellar CO radio line 
emission from a large sample of optically bright carbon stars, 
originally observed by Olofsson et al. (ApJS, 87, 267). Some 
new observational results are presented here. We determine 
some of the basic parameters that characterize circumstellar en- 
velopes (CSEs), e.g., the stellar mass loss rate, the gas expan- 
sion velocity, and the kinetic temperature structure of the gas. 
Assuming a spherically symmetric CSE with a smooth gas den- 
sity distribution, created by a continuous mass loss, which ex- 
pands with a constant velocity we are able to model reasonably 
well 61 of our 69 sample stars. The derived mass loss rates de- 
pend crucially on the assumptions in the circumstellar model, 
of which some can be constrained if enough observational data 
exist. Therefore, a reliable mass loss rate determination for 
an individual star requires, in addition to a detailed radiative 
transfer analysis, good observational constraints in the form of 
multi-line observations and radial brightness distributions. In 
our analysis we use the results of a model for the photodissoci- 
ation of circumstellar CO by Mamon et al. (1988). This leads 
to model fits to observed radial brightness profiles that are, in 
general, very good, but there are also a few cases with clear de- 
viations, which suggest departures from a simple r~ 2 density 
law. 

The derived mass loss rates span almost four orders of mag- 
nitude, from — 5xlO~ 9 MQyr -1 up to ~ 2x 10~ 5 M Q yr"\ 
with the median mass loss rate being 2.8x 10 _7 Mq yr~ 2 . We 
estimate that the mass loss rates are typically accurate to ^50% 
within the adopted circumstellar model. The physical condi- 
tions prevailing in the CSEs vary considerably over such a large 
range of mass loss rates. Among other things, it appears that the 
dust-to-gas mass ratio and/or the dust properties change with 
the mass loss rate. We find that the mass loss rate and the gas 
expansion velocity are well correlated, and that both of them 
clearly depend on the pulsational period and (with larger scat- 
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* Presented in this paper is observational data collected using the 
Swedish-ESO submillimetre telescope, La Silla, Chile, the 20 m tele- 
scope at Onsala Space Observatory, Chalmers Tekniska Hogskola, 
Sweden, and the NRAO 12 m telescope located at Kitt Peak, USA. 



ter) the stellar luminosity. Moreover, the mass loss rate corre- 
lates weakly with the stellar effective temperature, in the sense 
that the cooler stars tend to have higher mass loss rates, but 
there seems to be no correlation with the stellar C/O-ratio. We 
conclude that the mass loss rate increases with increased regu- 
lar pulsation and/or luminosity, and that the expansion velocity 
increases as an effect of increasing mass loss rate (for low mass 
loss rates) and luminosity. 

Five, of the remaining eight, sample stars have detached 
CSEs in the form of geometrically thin CO shells. The present 
mass loss rates and shell masses of these sources are estimated. 
Finally, in three cases we encounter problems using our model. 
For two of these sources there are indications of significant de- 
partures from overall spherical symmetry of the CSEs. 

Carbon stars on the AGB are probably important in return- 
ing processed gas to the ISM. We estimate that carbon stars 
of the type considered here annually return ~0.05Mq of gas 
to the Galaxy, but more extreme carbon stars may contribute 
an order of magnitude more. However, as for the total carbon 
budget of the Galaxy, carbon stars appear to be of only minor 
importance. 

Key words: Stars: AGB and post-AGB - Stars: carbon - cir- 
cumstellar matter - Stars: late-type - Radio lines: stars 



1. Introduction 

The carbon star phenomenon, i.e., the presence of stars for 
which the abundance of carbon exceeds that of oxygen in the 
atmospheres, occurs during the final evolutionary stage of low 
to intermediate mass stars (~1 — IOMq; Wallerstein & Knapp 
1998), and it is believed to be due to a dredge-up process driven 
by quasi-periodic He-shell flashes (Straniero et al. 1997). These 
stars, located on the asymptotic giant branch (AGB), lose copi- 
ous amounts of matter (^10~ 8 — 10 -4 Mq yr _1 ) in an intense 
stellar wind. In fact, it is the mass loss that determines the AGB 
lifetime and not the nuclear burning processes. Hence, a deter- 
mination of the mass loss characteristics, i.e., the dependence 
on time, mass, chemistry, metallicity, etc., is of utmost impor- 
tance for our understanding of late stellar evolution. 
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The mass loss creates a circumstellar envelope (CSE) of gas 
and dust around the star. The low temperature of the central star 
allows the formation of a wide variety of molecular species in 
its atmosphere, and the expanding gas is therefore mainly in 
molecular form. The chemistry in the CSE itself can be very 
rich, and it depends on the C/O-ratio, the thickness of the enve- 
lope, and the strength of the ambient ultraviolet radiation field. 
Presently, ^60 different molecular species have been identified 
in CSEs around AGB-stars (Olofsson 1997). 

The study of AGB-CSEs is important for the understanding 
of the late stages of stellar evolution as mentioned above. In ad- 
dition, an understanding of the atmospheric and circumstellar 
chemistry is required to determine the elemental abundances 
in these winds, which contribute to the chemical evolution of 
the interstellar medium (Busso et al. 1999). Finally, since AGB 
stars are the progenitors of planetary nebulae (PNe), the CSEs 
are very likely important ingredients in the formation of PNe 
(Kwok 1993). 

Observations of molecular millimetre-wave line emission 
have proven to be one of the best tools for studying the struc- 
ture, kinematics, and chemistry of the CSEs, and CO and OH 
observations have proven particularly useful for the determi- 



Table 1. Data on telescopes and receivers used. 



nation of accurate mass loss rates [see Olofsson (1996) and 
references therein]. 

While the overall picture of mass loss on the AGB is fairly 
well understood, some of the underlying physical principles are 
not. The prevailing theory, which needs to be thoroughly tested 
observationally, is that the mass loss occurs in a two stage pro- 
cess. First the pulsations of the star deposit energy in the atmo- 
sphere, leading to a considerably increased scale height, and 
sufficient matter at low enough temperatures for efficient dust 
formation. Radiation pressure on the dust grains accelerate the 
dust wind to its terminal expansion velocity in the second stage. 
The gas is momentum coupled, through collisions, to the dust, 
and it will be effectively dragged along. The mechanisms that 
cause the mass loss to vary on a wide range of time scales 
still need to be pinpointed. Eventually, the possibility of mass 
ejected in clumps has to be addressed, although in this paper 
we will consider only a smooth spherically symmetric wind. 

This study of an essentially complete sample of visually 
bright, relatively unobscured carbon stars will hopefully in- 
crease our knowledge of the modelling of circumstellar molec- 
ular radio line emission, the mass loss characteristics, the cir- 
cumstellar chemistry, the elemental abundances in the winds, 
and possibly the evolutionary status of these objects. It will also 
provide us with the tools required to determine the properties 
of higher mass loss rate objects, for which the central star may 
be completely obscured. 

In this paper the basic physical characteristics of the CSEs, 
e.g., the stellar mass loss rate, the gas kinetic temperature, 
and the gas expansion velocity, are determined using a ra- 
diative transfer analysis of the observed 12 CO (hereafter CO) 
millimetre-wave line emission. The CO molecule is particu- 
larly well suited for this purpose, since it is difficult to pho- 
todissociate and easy to excite through collisions, and thus is a 
very good tracer of the molecular gas density and temperature. 
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Furthermore, the molecular excitation is of a quasi-thermal na- 
ture which simplifies detailed modelling. 



2. Observations 

2.1. The sample 

We have selected a sample of 68 bright N- and J-type car- 
bon stars for which circumstellar CO emission was detected 
by Olofsson et al. (1993a). Their original sample consisted of 
carbon stars brighter than K—2 mag., in total 120 sources. The 
detection rate for all the carbon stars searched for circumstel- 
lar CO was as high as 72%. The original sample is thought to 
be close to complete to distances below 1 kpc [it is estimated 
that about a third of the carbon stars within this distance are 
missing, presumably the ones with thick and dusty CSEs, i.e., 
the high mass loss objects (Olofsson et al. 1993a)], and for dis- 
tances below 500 pc all sources (a total of 41) were detected 
in CO (note that the distances used in this paper are generally 
lower than those adopted in Olofsson et al. 1993a). These facts 
taken together led to the conclusion that the great majority of 
all N-type carbon stars are losing mass at a rate higher than 
about 10 _8 MQyr^ 1 . They do not, however, constitute a ho- 
mogeneous group in terms of circumstellar properties. Instead 
they show a wide range of gas expansion velocities and mass 
loss rates, presumably due to differences in mass, evolutionary 
stage, etc.. In addition, the reasonably well studied, high mass 
loss rate carbon star, LP And (also known as IRC +40540), is 
included in our analysis, since it provides a good test case for 
the radiative transfer model. 

The CO (J=l— >0 and J=2— *1) observations presented 



in Olofsson et al. (1993a) were obtained using the 20 m tele- 
scope at Onsala Space Observatory (OSO), Sweden, the 15 m 
Swedish-ESO sub-millimetre telescope (SEST) at La Silla, 
Chile, and the IRAM 30 m telescope at Pico Veleta, Spain, dur- 
ing the years 1986-1992. 

2.2. New observations 

We have supplemented the original data with new CO observa- 
tions of some of the sample stars. In particular, in the J=3— >2 
line using the SEST (August 1992). Also some J=1^0 and 
J=2^1 data were obtained at SEST (October 1998). A few 
sources were observed in the J=1^0 and J— 2— >1 lines with 
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Fig. 1. New observations of circumstellar 12 CO line emission. The line observed and the telescope used are shown in the upper 
right corner of each panel. 
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Table 2. Observational results (see text for details). The intensity scale is thought to be accurate to within ^20%. 
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1 Contamination by interstellar lines. 2 (3 was fixed. 



the NRAO 12 m telescope at Kitt Peak, USA (June 1998), and 
some J=1^0 data were obtained using OSO (1998 to 1999). 

A summary of telescope and receiver data [the receiver 
temperature T roc in single sideband mode (SSB), the full half 
power main beam width # m b, the main beam efficiency /7 m b, 
and the corrected main beam efficiency 77*J at the observational 
frequencies are given in Table [l]. 

At SEST, a dual channel SIS mixer receiver was used to 
simultaneously observe at 115 GHz (the J=l— >0 line) and 
230GHz (the J=2^1 line). At 345 GHz (the J=3^2 line), a 
single polarization SiS mixer receiver was used. As backends, 
two acousto-optical spectrometers (AOS) were used; one high 
resolution spectrometer (HRS) and one low resolution spec- 
trometer (LRS). The wideband (1 GHz) LRS had 1440 chan- 
nels separated by 0.7 MHz, whereas the narrow band (84 MHz) 
HRS used 2000 channels separated by 42 kHz. 

For the NRAO 12 m observations we used a dual polariza- 
tion SIS receiver. As backends two filterbanks, each with 256 
channels and a channel separation of 1 MHz, were used. The 
spectra obtained at each of the polarizations were added to in- 
crease the signal-to-noise ratio. In addition, for strong sources, 
a millimetre autocorrelator (MAC) was used in order to obtain 
a higher velocity resolution. The MAC was used in a configu- 
ration where the usable bandwith was 300 MHz and the resolu- 
tion was 98 kHz. 

At OSO, a single polarization SIS receiver was used for 
the observations. As backends two filterbanks, with bandwidths 



of 512 MHz (MUL A) and 64 MHz (MUL B), were used. 
The MUL A used 512 channels separated by 1 MHz, and the 
MUL B filterbank used 256 channels with a separation of 
250 kHz. 

The SEST and OSO observations were made in a dual 
beamswitch mode, where the source is alternately placed in the 
signal and the reference beam, using a beam throw of about 
12' at SEST and 11' at OSO. This method produces very flat 
baselines. The intensity scales are given in main beam bright- 
ness temperature, T m b. T m b=T£ / r] m ^, where T£ is the antenna 
temperature corrected for atmospheric attenuation using the 
chopper wheel method, and ry m b is the main beam efficiency, 
see Table [!} 

At the NRAO 12 m telescope the observations were carried 
out using a position switching mode, with the reference posi- 
tion located +10' in azimuth. This is the preferred observation 
mode for spectral line observations at the 12 m telescope. The 
raw spectra, which are stored in the scale, were converted 
using T m b= 2^/77^, where 77^ is the corrected main beam effi- 
ciency given in Table 0. The scale is related to T£ through 
Tfi=T£/T]f ss , where 77f ss is the forward scattering and spillover 
efficiency. 

Regular pointing checks were made on strong SiO masers 
(OSO, SEST) and strong continuum sources (NRAO). The 
pointing was usually consistent within ^3" for SEST and OSO 
and ^5" at NRAO. The uncertainties in the absolute inten- 
sity scales at the various telescopes are estimated to be about 
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±15—20%. However, due to the low efficiency of the SEST at 
345 GHz, and the narrow beam, we consider these data to be 
particularly uncertain in terms of calibration. 

The new observational results are summarized in Table [| 
and the spectra are shown in Fig. [j]. The line parameters, i.e., 
the main beam brightness temperature at the line centre (T m b), 
the line centre velocity and half the full line width (v e ), are 
obtained by fitting the following line profile to the data (Olofs- 
son et al. 1993a) 



T(v) = T mb 



1 



(3/2 



(1) 



where (3 is a parameter describing the shape of the line ((3=2 
represents a parabolic line shape, and (3<Q means a profile with 
horns at the extreme velocities). The integrated intensity (J m b) 
is obtained by integrating the emission between u*±u e . 

In addition, we have obtained publicly available data from 
the James Clerk Maxwell Telescope (JCMT) at Mauna Kea, 
Hawaii. The JCMT data are taken at face value. However, in the 
cases where there are more than one observation available, the 
derived line intensities are generally consistent within ±20%. 
The good agreement with corresponding SEST intensities fur- 
ther support the reliability of the JCMT public data. Further- 
more, interferometer observations of the CO( J=1^0) bright- 
ness distribution around some of our carbon stars have been 



obtained by Neri et al. (|1998|) using the IRAM Plateau de Bure 
interferometer (PdBI). The data are publicly available and have 
been used in this paper. 

3. Radiative transfer 

3.1. The Monte Carlo method 

In order to model the circumstellar molecular line emission and 
to derive the basic characteristics of the CSEs we have devel- 
oped a non-LTE radiati ve tran sfer code based on the M onte 
Carlo method [Bernes ( 1979 ); see also Schoier ( 2000 ), for 
details]. An accurate treatment of the molecular excitation is 
needed for a correct description of the radiative transfer in 
an expanding circumstellar envelope, e.g., the Sobolev ap- 
proximation has been shown to introduce significant errors in 
these circumstances (e.g., Schonberg 1985). The Monte Carlo 
method is well suited for the study of CSEs, since it is very 
flexible and close to the physics and yet simple in principle. For 
instance, one may include very complex geometries and veloc- 
ity fields without being forced away from the physics of the 
problem. Crosas & Menten 1997 have recently used the Monte 
Carlo method to model the circumstellar CO radio line emis- 
sion of the prominent carbon star IRC+10216 (from hereon 
called CW Leo). 

The aim of the radiative transfer is to obtain the steady-state 
level populations, of the molecule under study, using the statis- 
tical equilibrium equations (SE). In the Monte Carlo method 
information on the radiation field is obtained by simulating the 
line photons using a number of model photons, each represent- 
ing a large number of real photons from all transitions consid- 



ered. These model photons, emitted both locally in the gas as 
well as injected from the boundaries of the CSE, are followed 
through the CSE and the number of absorptions are calculated 
and stored. Photons are spontaneously emitted in the gas with 
complete angular and frequency redistribution (CRD), i.e., the 
local emission is assumed to be isotropic and the scattering are 
assumed to be incoherent. The weight of a model photon is 
continuously modified, to take the absorptions and stimulated 
emissions into account, as it travels through the CSE. When 
all model photons are absorbed in, or have escaped, the CSE 
the SE are solved and the whole process is then repeated until 
some criterion for convergence is fulfilled. Once the molecular 
excitation, i.e., the level populations, is obtained the radiative 
transfer equation can be solved exactly. 

The main drawback of the Monte Carlo method is its slow 
convergence (~v^itcr). This is, however, usually outweighed 
by its great adaptability. One of our motivations for choosing 
the Monte Carlo method was to be able to treat varying mass 
loss rates, departures from spherical symmetry, and the pos- 
sibility of a highly clumped medium. These complicating is- 
sues will be investigated in a future version of our Monte Carlo 
code. 

3.2. The standard model 

The numerical modelling of line emission from CSEs, where 
intricate interplays between physical and chemical processes 
take place, is a challenge. This applies, in particular, to the im- 
portant inner regions of the CSE model, where very few ob- 
servational constraints are available and our knowledge is very 
limited. In this analysis, we will neglect many of these com- 
plexities and assume a relatively simple, yet reasonably realis- 
tic, CSE model. By applying the model to a large sample we 
will derive fairly reliable mass loss rates, and we will also be 
able to pick out objects for which there appears to be clear de- 
viations from the simple picture. 

In what will be referred to as "the standard model" we as- 
sume a spherically symmetric CSE that expands at a constant 
velocity. In our data we see no evidence for gas acceleration 
in the CO emitting region. The CSE is divided into a number 
of discrete shells, each with its own set of physical parameters 
describing the state of the molecular gas. The density structure 
(as a function of distance, r, from the central star) can then be 
derived from the conservation of mass 



M 



PH 2 



n H2 m H2 



47rr 2 tv 



(2) 



where M is the hydrogen gas mass loss rate, and v c is the gas 
expansion velocity taken from the line profile fits of Olofsson 
et al. (1993a) or Table^. We assume that in the CSEs of interest 
here the hydrogen is in molecular form in the region probed by 
the CO emission (Glassgold & Huggins, 1983). The turbulent 
velocity is assumed to be equal to 0.5kms _1 throughout the 
entire CSE. 

The excitation of the CO molecules were calculated tak- 
ing into account 30 rotational levels in each of the ground and 
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first vibrational states. The radiative transition probabilities and 
energy levels are taken from Chandra et al. (1996), and the ro- 



tational collisional rate coefficients (CO-H2) are based on the 



results in Flower & Launay (1985). These are further extrapo 



lated for J> 1 1 and for temperatures higher than 250 K. Colli- 
sional transitions between vibrational levels are not important 
due to the low densities and their short radiative lifetimes. 

The kinetic gas temperature is calculated after each itera- 
tion in a self-consistent way, using the derived level popula- 
tions, by solving the energy balance equation (e.g., Goldreich 
& Scoville 1976), 



ar r nn 2 kv c 



(3) 



where 7 is the adiabatic index [assumed to be equal to 5/3, 
but the results are not very sensitive to its exact value, see also 
Groenewegen (|1994|) and Ryde et al. (|1999[)], k is the Boltz- 



mann constant, H is the total heating rate per unit volume, and 
C is the total cooling rate per unit volume. The first term on the 
right hand side is the cooling due to the adiabatic expansion of 
the gas. Additional cooling is provided by molecular line emis- 
sion from CO and H2. The molecular cooling due to CO is 
calculated from the derived level populations using the expres- 



sion of Crosas & Menten (1997 ). For the H2 cooling we use 
the approach by Groenewegen ( 1994 ). HCN could be an im- 
portant coolant in the inner par ts of t he envelope (Cernicharo 
et al. 1996), but Groenewegen ( |l 994| ) demonstrates that HCN 
cooling is only of minor importance. In the present version of 
the code HCN cooling is therefore not included (in a forth- 
coming paper, Schoier & Olofsson in prep., this issue will be 
addressed). 

The mechanism responsible for the observed mass loss is 
probably (at least for the stars of interest here) radiation pres- 
sure acting on small dust grains, which in turn are coupled to 
the gas (e.g., Hofner 1999). The radiation pressure on the dust 
grains will give them a drift velocity, «dr, relative to the gas 
(Gilman 1972; Goldreich & Scoville 1976; Kwan & Hill 1977) 



Lv c Q 
Mc 



1/2 



(4) 



where L is the luminosity, Q is the averaged momentum trans- 
fer efficiency, and c is the speed of light [see Sahai (1990) for a 
more elaborate treatment]. As a result of the dust-gas drift, ki- 
netic energy of the order of ^Tii^di w ^ ^ e transfered to the 
gas each time a particle collides with a dust grain. This is as- 
sumed to provide the dominating heating of the gas (Goldreich 
& Scoville 1976; Kwan & Hill 1977), 



ffdg = (nd<7d«dr)^PH 2 «dr! 



(5) 



where rid is the number density of dust grains, and Od the geo- 
metrical cross section of a dust grain. Introducing the dust-to- 
gas mass ratio if? this equation can be written 



Ha 



'H 2 "H 2 



1 'dr 



OdPd 1 + ^ 



(6) 



where ad and pd are the average size and density of a dust 
grain, respectively. 

When solving the energy balance equation free parameters 
describing the dust, i.e., via the dust-gas collisional heating, are 
introduced. These are highly uncertain, but affect the derived 
line intensities. Here we assume that the Q-parameter, i.e., the 
efficiency of momentum transfer, to be equal to 0.03 [for de- 
tails see Habing et al. (1994)], and define a new parameter that 
contains the other dust parameters, 



h : 



OOl 



2.0 gem 

Pd 



0.05 /im 



(7) 



where the normalized values are the ones we used to fit the CO 
line emission of CW Leo using our model, i.e., h=l for this 
object. 

Additional heating is provided by the photoelectric effect 
(Huggins et al. 1988) 



(8) 



where /c pc is a constant that depends on the properties of 



the dust grains. Following Huggins et al. (1988) we adopt 



£; pe =lxl0 ergs . This heating mechanism is important 
in the cool, tenuous, outer parts of CSEs around high mass loss 
rate stars. 

The spatial extent of the molecular envelope is generally 
an important input parameter, and the derived mass loss rate 
will depend on this. The radial distribution of CO in the CSE 
was es timated using the modelling presented in Mamon et al. 
( 1988 ). It includes photodissociation, taking into account the 
effects of dust-, self- and ^-shielding, and chemical exchange 



reactions. Mamon et al. (|1988j) find that the radial abundance 
distribution of CO with respect to H2, f(r), can be described 
by 



f(r) = fo exp 



-ln2 



(9) 



where /o is the initial (photospheric) abundance, r p is the 
photodissociation radius (where the abundance has dropped to 
fo/2), and a is a parameter describing the rate at which the 
abundance declines. Both r p and a depend on the mass loss 
rate (M), the expansion velocity (v c ), and fo. In our mo dellin g 
we use analytical fits to these results from Stanek et al. ( 1995 ) 



5.4 x 10 16 



M 
HP 



-0.55 



fo 



8 x 10" 4 



and Kwan & Webster ( 1993 ) 



a = 2.79 



M 
1(P 



0.09 



0.09 



0.55 



(10) 



(11) 



where the units of M and v c are MQyr 1 and km s 1 , respec- 
tively. In what follows we will assume /o=l x 10~ 3 . This is an 
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average of the /o:s estimated by Olofsson et al. ( 1993b| ) for 
this sample of optically bright carbon stars. The scatter in the 
estimated /o :s 1S about 40%. We note here that a somewhat 
more sophisticated photodissociation model was developed by 



Doty & Leung (1998) in which they include scattering by dust. 



For the CSE of CW Leo they derive a CO envelope size which 
is 30% smaller than what is obtained using the Mamon et al. 
model. It is however not clear by how much the envelope size 
changes for an object with a significantly lower mass loss rate, 
and hence we use the results of Mamon et al. here. 

In our models we include both a central source of radia- 
tion and the cosmic microwave background radiation at 2.7 K. 
The central radiation emanates from the star, and was estimated 
from a fit to its spectral energy distribution (SED), usually by 
assuming two blackbodies. This method is described in Ker- 
schbaum fit to the SED gives the two blackbody 

temperatures, T* and Td, and the relative luminosities of the 
two blackbodies, L^IL*. Any dust present around the star will 
absorb parts of the stellar light and re-emit it at longer wave- 
lenghts. Thus, of the two blackbodies used, one represents the 
stellar contribution and one represents the dust. However, it 
should be noted that, e.g., the stellar temperature T» derived 
in this manner is generally lower than the true temperature 
of the star (Kerschbaum 1999). In any case, the two black- 
bodies used give a good description of the radiation that the 
CSE is subjected to. The temperatures and luminosities used 
in the modelling are presented in Table ^. Only in the cases 
where Ld/£*>0.1 the dust component was retained. The inner 
boundary of the CSE, r\, was set to reside outside the radius 
of the central blackbody (s), but never lower than lxl0 14 cm 
(~3Rq), i.e., generally beyond both the sonic point and dust 
condensation radius. 

The distances, presented in Table ||, were estimated using 
one of the following methods: the observed Hipparcos paral- 
lax, a period-luminosity relation (Groenewegen & Whitelock 
1996), or an assumed bolometric luminosity. In the two for- 
mer cases the luminosities were estimated using apparent bolo- 
metric magnitudes and the distances. In the rare cases when 
both the distance obtained from Hipparcos and the period- 
luminosity relation result in extremely high or low luminosi- 
ties, a value of 4000 Lq was adopted as the luminosity and the 
distance was estimated using this value. For the irregular vari- 
ables with no Hipparcos data we also assumed a luminosity 
of 4000 Lq. For a statistical study of a large sample of stars 
these distance estimates are adequate, although the distance es- 
timate for an individual star has a large uncertainty. There is no 
apparent systematic difference between the distances derived 
from the period-luminosity relation and the Hipparcos paral- 
laxes, Fig. ||. The distances in Olofsson et al. (1993a) were es- 
timated from an adopted absolute K magnitude. The distances 
presented here are systematically lower by, on the average, a 
factor of 1 .4 than those used in Olofsson et al. 

This was a summary of the assumptions used in "the stan- 
dard model". In what follows, tests will be made to see how 
sensitive the model is to the various assumptions. 



o 
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Fig. 2. Comparison between the distances derived from the 
period-luminosity relation (Dp_l) and the Hipparcos paral- 
laxes (-Dhip)- The irregular variables where the luminosity was 
assumed to equal 4000 Lq are indicated by open circles. The 
dashed line shows the 1 : 1 correlation. 



4. Model results 
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Fig. 3. The ^.-parameter derived from the radiative transfer 
analysis plotted against the adopted luminosity of the central 
star. Note that the data points at L=4000Lq represent two 
stars each. 
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Table 3. Input parameters and CO modelling results using the standard model (see text for details) 



Source 


Var. type 


P 


D 


L 


T* 






[days] 


[pc] 


[L©] 


[K] 


v a /\nci 


CDn 

oKa 


^AQ 

jOy 


jOU 


3 JUU 


zuuu 


1 1 V l_d» 


A/T 

1V1 


S97 

JZ 1 


Q70 2 

y /u 




loUU 


TP C _i_AnnA 1 
1KC+0UU4 1 






/u 


4UUU 


1 ^nn 


Z, r SC 


CD], 


1 AA 
144 


SIX) 


^ 1 nn 

J 1UU 


ZOUU 


P T7nr 


A/T 
1V1 


J07 


Olv^ 


JoUU 


1 JUU 


1 W 1 101 


QPVi 


1 S8 

1 Jo 


Ann 1 


A7nn 
o /uu 


o Ann 

ZOUU 


V38A Dor 

v jo4 r er 


A/T 




jOU 


s 1 nn 


1 onn 
ryuu 


VlAfi Par 

v "400 r er 


CD 

jK 






4UUU 


oonn 
zzuu 


o i Lam 




^nn 

JUU 


^n 2 


AAOO 
^-^fUU 


1 Ann 

10UU 


r i 1 Ton 
i i iau 


QPVi 


1 fn 
10 / 


JOU 


9Ann 

Z4UU 


OAnn 

Z4UU 


P T on 

k Lep 


A/T 


AA7 

40 / 


9^n! 

LJ\) 


Annn 

4UUU 


1 onn 
ryuu 


w un 


CD], 


919 

Z1Z 


99H 1 
ZZU 


9 Ann 
zouu 


oonn 
zzuu 


Q A iii- 

o /\ur 


en 


OVU 


(5ZU 


QQflA 


1 Ann 
louu 


W Pi,- 


T h 






Annn 


?nnn 
zuuu 


V Tan 

i iau 




9A9 
Z4Z 


^ 1 n 2 


j juu 


OAnn 

Z^fUU 


1 U v J c FT) 




9^n 

ZjU 


jjU 


^Ann 

J4UU 


9Ann 

Z^fUU 


P.T Ovi 


T l-i 






Annn 


97nn 
z / uu 


TTTT Anr 




9X1 
Zj4 


9AH 2 

zou 


Aonn 
oyuu 


osnn 

ZJUU 


MP Pun 


T Vi 

LD 




A9H 1 
"4-ZU 


^snn 
j juu 


97nn 
z / uu 


ll ivion 


A/T 


AQ7 

4y / 


77H 2 


7^nn 
/ juu 


onnn 
zuuu 


PV A/Tnn 


ODn 

oKil 


ASA 
4JU 


ASA 1 


^nnn 


9nnn 
zuuu 


w LJvia 


T l-i 






/innn 
4UUU 


osnn 
ZjUU 


P Vnl 
K VOl 


A/T 

1V1 


ASA 


7^n 2 


Asnn 

OoUU 


i ^nn 

1 JUU 


A 1_T1C 




1 QS 

lyD 


9sn 2 

ZoU 


9snn 

ZoUU 


oonn 
zzuu 


lw Leo 


A/T 
JV1 


A'in 
OJU 


1 9fl 2 
1ZU 


QAnn 
V0UU 




i ny a 




jUj 


jjU 


A9nn 


9 Ann 

ZOUU 


V Vol 

a vei 


cp 


1 AH 
14U 


^ 1 n 2 


9snn 

ZoUU 


oonn 
zzuu 


o£ Car 


cnu 


1 9A 
1Z0 




/innn 
4UUU 


o/i nn 
Z4UU 


PW T A/Ti 




AAn 

04U 


AAO 2 


Q7nn 
y / UU 


1 ^nn 

1 jUU 


Y7 Vol 

al vei 








4UUU 


1 Qnn 
iyuu 


ca riya 


A/T 


/1/19 
44Z 


you 


AAnn 
OoUU 


1 Qnn 


PPC 97Q9 






A so 3 


AC\(\C\ 
H-UUU 


9nnn 
zuuu 


u riya 


QPVi 
oKD 


Asn 

4jU 


1 AO 1 
10U 


9^nn 

ZjUU 


9Ann 

Z4-UU 


VV T TA/Ta 

v i uivia 


T U 




jjV 


4UUU 


97nn 
z / uu 


vir 




^AA 


SAO 2 


SAnn 


1 snn 

loUU 


V fVn 


QPVi 
oKD 


1 S7 


99H 1 
ZZU 


44UU 


9?nn 
zzuu 


PV Plro 

Ki ura 


cnu. 
oKD. 


9nn 
zuu 


/ion 1 
4VU 


1 nnnn 
1UUUU 


osnn 

ZjUU 


T4r> 1 9 1 asr 
1 1 1 J iziDJo 






JOU 


"4UUU 


Z4UU 


l in 1 9A9AS 
1 1 l_y rZH-ZOo 






JVU 




9^nn 

ZJUU 


V Tr A 
A 1 1 A 


T Vi 




aao 3 

40U 


4UUU 


99nn 
zzuu 


V r r t> 

v ltd 


Ayr 




A^n 2 

OjU 


s^nn 

J jUU 


1 Qnn 
1 vuu 


i w vjpn 


QPVi 
oKD 


1 ss 

loJ 


9SH 1 
ZoU 


97nn 
z / uu 


onnn 
zuuu 


i ura 


A/T 


A99 
4ZZ 


A10 2 

01U 


A^nn 

OjUU 


1 Ann 

10UU 


T T \rr 


T h> 






AC\(\C\ 
tuuu 


9nnn 
zuuu 


V A /,| 


QPVi 
oKD 




^70 1 
J /U 


A^nn 
ojuu 


1 1 nn 

Z1UU 


V1Q4? S(rr 


Lb 




430 3 


4000 


2600 


UXDra 


SRa: 


168 


310 3 


4000 


2600 


AQSgr 


SRb 


200 


420 3 


4000 


2700 


RT Cap 


SRb 


393 


450 2 


5900 


2200 


UCyg 


M 


463 


710 2 


6900 


2200 


VCyg 


M 


421 


370 2 


6200 


1500 


RV Aqr 


M 


454 


670 2 


6800 


1300 


TInd 


SRb 


320 


570 1 


8800 


2600 


YPav 


SRb 


233 


360 2 


4100 


2400 


V1426 Cyg 


M 


470 


780 2 


7100 


1200 


SCep 


M 


487 


340 2 


7300 


1400 







M 






h 


[K] 




[MeyA] 


[kmr 1 ] 


[cm] 








4.0x 10~ 8 


11.5 


1.4x 


10 16 


0.2 4 


900 


0.54 


9.0xl0 -7 


18.5 


6.1x 


10 16 


0.5 4 


500 


0.13 


3.0x 10~ 6 


28.0 


l.Ox 


10 17 


0.2 4 






2.5x 10 -8 


3.5 


1.4x 


LO 16 


0.2 4 


700 


0.52 


1.3x 10 -6 


16.5 


8.0x 


10 16 


0.7 






9.0x 10 _s 


5.5 


2.5 x 


LO 16 


0.5 


900 


10 


3.5xl0~ 6 


15.0 


1.4x 


10 17 


0.4 






1.3x 10 -7 


9.0 


2.6x 


LO 16 


0.2 4 






9.0x 10 -8 


9.0 


2.1x 


LO 16 


0.2 4 






7.0x 10~ 8 


5.0 


2.2x 


LO 16 


0.2 4 


800 


0.14 


7.0x 10 -7 


18.0 


5.3x 


10 16 


0.2 






7.0x 10 -8 


11.0 


1.8x 


LO 16 


0.3 


700 


0.32 


2.2x 10 -6 


25.5 


8.9x 


10 16 


0.5 4 






3.0xl0~ 7 


16.0 


3.5x 


to 16 


0.4 






4.0x 10 -7 


11.0 


4.5 x 


LO 16 


0.1 






2.5xl0 -7 


11.5 


3.4x 


10 16 


0.2 4 






l.lxl0~ 7 


9.0 


2.4 x 


LO 16 


0.2 4 






3.5xl0 -7 


11.0 


4.2x 


LO 16 


0.1 






6.5 xlO -8 


9.5 


1.8x 


L0 16 


0.2 4 


1100 


1.4 


2.2xl0~ 6 


25.0 


8.9x 


LO 16 


0.5 4 






3.5xl0 -7 


11.0 


4.2x 


LO 16 


0.2 4 






3.0x 10 -7 


10.5 


3.9x 


LO 16 


0.2 4 


800 


0.7 


1.8xl0 -6 


18.0 


9.0x 


LO 16 


1.0 






l.Ox 10 -7 


7.0 


2.4 x 


LO 16 


0.1 


510 




1.5xl0 -5 


14.5 


3.7x 


L0 17 


1.0 






1.9xl0~ 7 


9.0 


3.2x 


LO 16 


0.2 4 






1.8x 10 -7 


10.0 


3.0x 


LO 16 


0.2 4 






2.0xl0 -6 


14.0 


l.lx 


L0 17 


0.2 


510 


6.7 


6.0xl0 -6 


17.0 


1.9x 


L0 17 


1.4 






6.0x 10~ 7 


14.0 


5.2x 


LO 16 


0.2 4 






9.0x 10 -7 


12.0 


7.0x 


LO 16 


0.5 4 






6.5xl0 -7 


17.0 


5.2x 


LO 16 


0.2 4 






1.4xl0 -7 


7.0 


2.9x 


LO 16 


0.2 4 






7. Ox 10 -8 


6.0 


2.1x 


10 16 


0.2 4 






2.0x 10 -7 


12.5 


3.0x 


LO 16 


0.2 4 






1.5x 10~ 7 


8.5 


2.9x 


LO 16 


0.2 5 






3.0x 10 -7 


10.0 


4.0x 


LO 16 


1.2 5 






l.OxlO -7 


6.5 


2.5 x 


L0 16 


0.2 4 






l.Ox 10~ 7 


11.0 


2.2x 


L0 16 


0.2 4 






1.3x 10~ 7 


8.0 


2.7 x 


LO 16 


0.4 


740 


0.15 


6.0xl0 -7 


7.5 


6.8x 


LO 16 


0.2 






5.0xl0 -8 


7.5 


1.6x 


LO 16 


0.2 4 


650 


0.28 


1.2x 10~ 6 


13.5 


8.0x 


LO 16 


0.5 4 






7.0xl0 -8 


11.5 


1.8x 


10 16 


0.5 5 






3.0x 10 -7 


8.5 


4.2x 


10 16 


0.2 






1.6x 10 -7 


in n 


2.8x 


10 16 


n ? 4 






1.6xl0 -7 


4.0 


4.0x 


LO 16 


0.2 4 






2.5 xlO -7 


10.0 


3.6x 


LO 16 


0.2 4 






l.OxlO -7 


8.0 


2.3 x 


LO 16 


0.2 4 






9.0xl0 -7 


13.0 


6.8x 


LO 16 


0.5 4 


860 


0.84 


1.2xl0 -6 


11.5 


8.5x 


LO 16 


1.2 


620 


0.46 


2.5xl0~ 6 


16.0 


l.lx 


L0 17 


0.5 






9.0xl0 -8 


6.0 


2.4 x 


LO 16 


0.5 4 






1.6xl0 -7 


8.0 


3.0x 


LO 16 


0.2 4 


580 


0.5 


l.OxlO -5 


14.0 


2.9x 


L0 17 


0.3 






1.5xl0~ 6 


22.0 


7.5x 


10 16 


0.4 
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Table 3. continued. 



Source Van type P D L T, T d 6 L d /Ll M v c r p h 

[days] [pc] [L Q ] [K] [K] [Moyr" 1 ] [kms^ 1 ] [cm] 



V460 Cyg 


SRb 


180 


230 2 


2600 


2800 






1.8x10" 


7 


10.0 


3.0xl0 16 


0.2 4 


RVCyg 


SRb 


263 


310 2 


3400 


2100 






4.5x10" 


7 


13.5 


4.5 xlO 16 


0.2 4 


PQ Cep 


M 




390 2 


4000 


1700 


840 


0.13 


1.4x10" 


6 


19.5 


7.6xl0 16 


0.2 4 


LP And 


M 


620 


630 2 


9400 


1100 


610 


6.6 


1.5x10" 


5 


14.0 


4.0xl0 17 


0.7 


WZCas 


SRb 


186 


290 2 


2700 


2500 






6.5x10" 


9 


2.5 


7.4xl0 15 


0.2 4 



1 Hipparcos; 2 Period-luminosity relation (Groenewegen & Whitelock 1996); 3 L=4000Lq assumed; 4 Assumed /i-parameter (see text for 
details); 5 J-type star (only including 12 CO cooling); 6 In the cases when L^/ L„<Q.l the dust blackbody component was not retained and a 
single blackbody representing the stellar component adequately describes the SED. In the case of CW Leo a single dust blackbody is used to 
represent the SED. 



4.1. The fitting strategy 

With the assumptions made in the standard model there re- 
mains two principal free parameters, the mass loss rate (M) 
and the /i-parameter. These two parameters were allowed to 
vary simultaneously, M in steps of ~10% and h in somewhat 
larger steps, until the model with the smallest deviations from 
the observed intensities was found. The quality of a particu- 
lar model with respect to the observational constraints can be 
quantified using the chi-square statistic defined as 



X 



N 

E 



[4 



^obs,i)] 



(12) 



where / is the total integrated line intensity, cr, is the uncer- 
tainty in observation i, and the summation is done over all in- 
dependent observations N. The errors in the observed inten- 
sities are generally dominated by the calibation uncertainty of 
^20%. In some cases the line profiles were used to discrimi- 
nate between models. For each new M assigned, a consistent 
envelope si ze (r p ) was calculated based upon the results of Ma- 
mon et al. (11988b. 



It is found that the intensity ratios between the different 
lines are sensitive to various parameters, and hence can be 



used to constrain them (Sec. 4.3). Therefore, in order to ac 



curately determine h, and ultimately M, multi-transition ob- 
servations are needed. Observed radial brightness profiles pro- 
vide additional important information. For a limited number of 
the sources in our sample there are enough observational con- 
straints that a reliable /i-parameter may be estimated. These 
sources were first modelled, see Table || 

4.2. The h-parameter 

The derived ^.-parameters, as a function of the adopted lumi- 
nosity of the central source, are shown in Fig. |[ A division 
between low luminosity objects and those with higher lumi- 
nosities is evident. For luminosities lower than 6000 Lq the 
median value of the estimated /i-parameter is 0.2 (based on 9 
sources), while for the higher luminosities it is 0.5 (12 sources) 
(the J-type stars, i.e., stars with a high content of 13 CO, were 
not included in these estimates since 13 CO line cooling is not 
treated in the standard model). The large scatter for the high 



luminosity sources is unfortunate, since these are generally the 
sources where the derived mass loss rates are particularly sen- 
sitive to the adopted value of h. For the majority of objects, 
where a reliable estimate of the ^-parameter is not possible, we 
have adopted a value of h~0.2 for the low luminosity objects 
(below 6000 Lq) and h—Q.5 for the high luminosity objects 
(above 6000 L ). 

The /i-parameter is related to the dust-to-gas mass ratio, if? 
(see Eq. fy. If we assume that the properties of a dust grain are 
equal to those used to fit the line emission from CW Leo, and 
that they are the same for all stars, an ^.-parameter of 0.2 means 
that \]>=2xl0~ 3 . Groenewegen et al. (|1998bh have modelled 



the dust emission towards a large sample of carbon stars and 
find a value of 'J/ which is fairly constant at 2.5 xlO -3 for 
stars with luminosities up to about 7900 Lq, and which in- 
creases drastically with luminosity up to 0.01 and higher. This 
is in good agreement with our results and lends further sup- 
port to our distinction between low- and high-luminosity stars 
as rega rds the value of the /i-parameter. Furthermore, Hiriart & 
Kwan (2000) derive, for 17 of our sample stars, on the aver- 
age ^~7xl0"" 4 . The main reason for this discrepancy is the 
difference in the adopted dust parameters. In fact, it can be 
shown that for h—0.2 we get essentially the same heating term 
(cf. Eq. ||) in the energy balance equation as Hiriart & Kwan 
fl200Cp , 



4.3. The mass loss rates 

We have managed to model reasonably well 61 of the 69 sam- 
ple sources using our standard model and fitting strategy de- 
fined above. The reduced \ 2 for these models is estimated from 



v 2 

2 A min 

Arcd 



N- 



(13) 



where Xmin i s obtained from Eq. |12J for the best fit model, 
and p is the number of adjustable parameters (one or two) in 
the model. Typically, x 2 od ~l— 3 assuming a calibration uncer- 
tainty of 20%. The derived mass loss rates and /i-parameters 
are presented in Table ^|, where also the velocities by which the 
CSEs expand, v c , and their spatial extents, r p , are given. 

The derived mass loss rates presented in Table || span 
more than three orders of magnitude, 6.5 x 10 _9 Mq yr _1 to 
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Table 4. CO modelling results compared to single dish observations. 



Source 


Tel. 


Trans. 


Zobs 

[Kkmr 1 ] 


•^raod 

[Kkms -1 ] 


Ref. 


Source 


Tel. 


Trans. 


-^obs 

[Kkms -1 ] 


-^mod 

[Kkms -1 ] 


Ref. 


VXAnd 


oso 


1-0 


<0.8 


0.1 


2 


UUAur 


OSO 


1-0 


7.9 


8.2 


2 




IRAM 


1-0 


<1.4 


0.2 


2 




IRAM 


1-0 


18.8 


16.9 


2 




IRAM 


2-1 


2.4 


2.4 


2 




JCMT 


2-1 


16.8 


14.8 


4 


HV Cas 


OSO 


1-0 


<2.1 


1.9 


2 




IRAM 


2-1 


39.0 


46.2 


2 




IRAM 


1-0 


5.6 


4.1 


2 


NP Pup 


SEST 


1-0 


<0.5 


0.3 


2 




IRAM 


2-1 


11.9 


14.5 


2 




SEST 


2-1 


1.2 


1.2 


2 


IRC+60041 


IRAM 


1-0 


IS 


4.7 


2 


CL Mon 


SEST 


1-0 


3.2 


3.3 


2 




IRAM 


2-1 


17.3 


17.3 


2 




SEST 


2-1 


10.4 


9.5 


2 


ZPsc 


OSO 


1-0 


1.0 


0.7 


2 




IRAM 


2-1 


34.3 


38.1 


2 




IRAM 


1-0 


2.1 


1.6 


2 


RYMon 


IRAM 


1-0 


6.8 


5.4 


2 




IRAM 


2-1 


4.6 


6.0 


2 




IRAM 


2-1 


11.3 


14.5 


2 


R For 


SEST 


1-0 


4.8 


4.9 


2 


WCMa 


SEST 


1-0 


1.3 


1.3 


2 




SEST 


2-1 


11.0 


14.5 


2 




SEST 


2-1 


4.6 


3.9 


2 




JCMT 


2-1 


21.0 


17.3 


4 




IRAM 


2-1 


11.7 


16.0 


2 




SEST 


3-2 


18.0 


19.9 


1 


R Vol 


SEST 


1-0 


4.9 


5.1 


2 




JCMT 


3-2 


27.4 


25.6 


4 




SEST 


2-1 


15.3 


14.7 


2 


TWHor 


SEST 


1-0 


1.0 


0.9 


2 




SEST 


3-2 


21.0 


22.2 


1 




SEST 


2-1 


3.3 


3.8 


2 


XCnc 


NRAO 


1-0 


0.7 


0.7 


1 




SEST 


3-2 


7.2 


6.1 


1 




SEST 


1-0 


1.5 


1.1 


2 


V384 Per 


NRAO 


1-0 


7.8 


10.1 


1 




OSO 


1-0 


1.7 


2.0 


2 




OSO 


1-0 


25.5 


25.1 


2 




SEST 


2-1 


3.2 


2.9 


2 




IRAM 


1-0 


52.7 


48.4 


2 




IRAM 


2-1 


11.1 


11.8 


2 




JCMT 


2-1 


37.3 


35.2 


4 


CWLeo 


NRAO 


1-0 


170.8 


257.1 


1 




IRAM 


2-1 


83.6 


102.0 


2 




SEST 


1-0 


288.1 


305.3 


2 




JCMT 


3-2 


44.5 


42.5 


4 




OSO 


1-0 


422.0 


391.7 


2 


V466 Per 


IRAM 


1-0 


3.6 


2.5 


2 




SEST 


2-1 


487.3 


580.5 


2 




IRAM 


2-1 


6.2 


8.6 


2 




JCMT 


2-1 


731.8 


632.2 


4 


ST Cam 


OSO 


1-0 


2.0 


2.0 


2 




IRAM 


2-1 


1057.7 


1116.1 


2 




JCMT 


2-1 


3.4 


4.4 


4 




SEST 


3-2 


854.8 


774.0 


1 




IRAM 


2-1 


15.2 


14.8 


2 




JCMT 


3-2 


1066.3 


883.7 


4 


TTTau 


IRAM 


1-0 


2.6 


2.7 


2 




JCMT 


4-3 


1227.8 


1052.0 


4 




IRAM 


2-1 


4.0 


8.1 


2 


YHya 


SEST 


1-0 


1.6 


1.6 


2 


RLep 


SEST 


1-0 


6.2 


6.7 


2 




IRAM 


2-1 


18.7 


19.4 


2 




IRAM 


1-0 


31.7 


24.3 


2 


X Vel 


SEST 


1-0 


1.4 


1.5 


2 




SEST 


2-1 


18.1 


17.4 


2 




SEST 


2-1 


4.3 


4.4 


2 




IRAM 


2-1 


56.4 


64.2 


2 


SZ Car 


SEST 


1-0 


2.8 


2.7 


2 




SEST 


3-2 


12.5 


21.6 


2 




SEST 


2-1 


7.8 


5.9 


2 




JCMT 


3-2 


28.1 


27.9 


4 




SEST 


3-2 


5.9 


6.2 


1 


WOri 


SEST 


1-0 


1.2 


1.3 


2 


RWLMi 


NRAO 


1-0 


36.4 


37.0 


1 




OSO 


1-0 


2.7 


2.4 


2 




SEST 


1-0 


54.6 


51.9 


2 




IRAM 


1-0 


5.0 


5.1 


2 




OSO 


1-0 


87.0 


83.8 


2 




SEST 


2-1 


4.9 


4.3 


2 




SEST 


2-1 


105.4 


128.4 


2 




IRAM 


2-1 


19.0 


16.9 


2 




JCMT 


2-1 


163.2 


147.4 


4 




SEST 


3-2 


4.8 


5.4 


1 




JCMT 


3-2 


245.3 


208.6 


4 


S Aur 


OSO 


1-0 


6.7 


5.5 


2 




JCMT 


4-3 


243.1 


224.3 


4 




IRAM 


1-0 


10.9 


11.9 


2 


XZ Vel 


SEST 


1-0 


1.8 


1.4 


2 




IRAM 


2-1 


15.8 


36.8 


2 




SEST 


2-1 


3.5 


4.1 


2 


WPic 


SEST 


1-0 


0.9 


0.9 


2 


CZ Hya 


SEST 


1-0 


1.5 


1.5 


2 




SEST 


2-1 


2.6 


3.1 


2 




SEST 


2-1 


4.2 


4.5 


2 




SEST 


3-2 


6.1 


4.5 


1 


CCS 2792 


SEST 


1-0 


1.9 


1.8 


2 


YTau 


SEST 


1-0 


3.2 


3.0 


2 




SEST 


2-1 


4.1 


4.7 


2 




OSO 


1-0 


6.3 


5.3 


2 


UHya 


SEST 


1-0 


5.4 


5.4 


2 




IRAM 


1-0 


15.5 


11.0 


2 




SEST 


2-1 


13.8 


14.1 


1 




SEST 


2-1 


5.8 


6.1 


2 




JCMT 


2-1 


20.2 


16.6 


4 




IRAM 


2-1 


19.4 


24.1 


2 




IRAM 


2-1 


48.8 


48.7 


4 


TU Gem 


OSO 


1-0 


3.4 


3.4 


2 


VYUMa 


OSO 


1-0 


1.4 


1.4 


2 




IRAM 


1-0 


6.5 


21.1 


2 




IRAM 


2-1 


4.3 


11.0 


2 


BL Ori 


OSO 


1-0 


IS 


1.2 


2 


SS Vir 


OSO 


1-0 


<1.2 


1.1 


2 




IRAM 


1-0 


IS 


2.6 


2 




IRAM 


1-0 


7.9 


8.5 


2 




IRAM 


2-1 


9.8 


9.7 


2 
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Table 4. continued. 



Source 


Tel. 


Trans. 


-^obs ^mod 


Ref. 


Source 


Tel. 


Trans. 


-^obs 


-^mod 


Ref. 








[Kkms -1 ] [Kkms -1 ] 










[Kkms -1 ] 


[Kkms -1 ] 




YCVn 


OSO 


1-0 


4.5 


5.4 


2 


VCyg 


OSO 


1-0 


27.5 


27.2 


2 




IRAM 


1-0 


10.3 


11.2 


2 




NRAO 


2-1 


50.1 


36.5 


1 




JCMT 


2-1 


11.9 


10.9 


4 




JCMT 


2-1 


35.9 


56.2 


4 




IRAM 


2-1 


22.7 


34.2 


2 




JCMT 


3-2 


88.9 


86.1 


4 




JCMT 


3-2 


20.9 


15.6 


4 




JCMT 


4-3 


123.4 


95.2 


4 


RYDra 


OSO 


1-0 


2.4 


3.0 


2 


RV Aqr 


SEST 


1-0 


7.5 


6.9 


2 




IRAM 


2-1 


21.6 


28.3 


2 




SEST 


2-1 


18.1 


17.2 


2 




JCMT 


3-2 


18.9 


15.8 


4 




SEST 


3-2 


18.6 


20.5 


1 


HD 121658 


SEST 


1-0 


1.5 


1.4 


2 


TInd 


SEST 


1-0 


0.5 


0.4 


2 


HD 124268 


SEST 


1-0 


<0.7 


0.6 


2 




SEST 


2-1 


1.5 


1.9 


2 




SEST 


2-1 


2.1 


1.9 


2 


YPav 


SEST 


1-0 


1.6 


1.4 


2 


XTrA 


SEST 


1-0 


2.5 


2.7 


2 




SEST 


2-1 


3.8 


4.3 


2 




SEST 


2-1 


11.8 


10.2 


2 


V1426 Cyg 


OSO 


1-0 


17.3 


17.0 


2 




SEST 


3-2 


15.3 


15.1 


1 




SEST 


2-1 


8.7 


14.3 


2 


VCrB 


OSO 


1-0 


2.7 


3.9 


2 




SEST 


3-2 


17.8 


17.1 


1 




IRAM 


1-0 


11.4 


8.4 


3 


SCep 


OSO 


1-0 


16.3 


19.4 


2 




IRAM 


2-1 


18.4 


18.8 


2 




IRAM 


1-0 


48.5 


38.1 


3 


TWOph 


SEST 


1-0 


1.1 


0.9 


2 




JCMT 


2-1 


49.3 


37.3 


4 




IRAM 


2-1 


8.7 


10.2 


2 




IRAM 


2-1 


76.6 


109.6 


2 


TDra 


OSO 


1-0 


8.8 


8.4 


2 




JCMT 


3-2 


55.7 


55.4 


4 




IRAM 


1-0 


28.6 


17.7 


3 




JCMT 


4-3 


65.7 


60.0 


4 




IRAM 


2-1 


43.7 


48.7 


2 


V460 Cyg 


OSO 


1-0 


4.7 


4.8 


2 


TLyr 


OSO 


1-0 


0.7 


1.0 


2 




IRAM 


2-1 


8.5 


28.4 


2 




IRAM 


2-1 


5.6 


9.6 


2 


RVCyg 


OSO 


1-0 


4.8 


6.1 


2 




JCMT 


3-2 


6.4 


4.3 


4 




IRAM 


1-0 


16.7 


12.9 


2 


V Aql 


SEST 


1-0 


2.8 


2.6 


2 




IRAM 


2-1 


14.4 


33.5 


2 




OSO 


1-0 


3.2 


4.8 


2 


PQCep 


OSO 


1-0 


6.4 


6.7 


2 




SEST 


2-1 


8.1 


7.0 


2 




NRAO 


2-1 


8.7 


7.4 


1 




JCMT 


2-1 


9.0 


8.3 


4 


LP And 


OSO 


1-0 


57.0 


58.8 


1 




JCMT 


3-2 


11.2 


10.6 


4 




IRAM 


1-0 


92.7 


97.4 


3 


VI 942 Sgr 


SEST 


1-0 


0.8 


0.8 


2 




NRAO 


2-1 


60.7 


46.5 


1 


UXDra 


OSO 


1-0 


2.1 


1.7 


2 




JCMT 


2-1 


90.6 


68.5 


4 




IRAM 


2-1 


6.8 


9.0 


2 




IRAM 


2-1 


155.8 


166.8 


3 


AQSgr 


SEST 


1-0 


1.4 


1.3 


2 




JCMT 


3-2 


88.0 


79.1 


4 




SEST 


2-1 


4.1 


4.0 


2 




JCMT 


4-3 


73.3 


78.6 


4 


RT Cap 


SEST 


1-0 


0.6 


0.6 


2 


WZCas 


OSO 


1-0 


<1.1 


0.1 


2 




SEST 


2-1 


2.0 


2.1 


2 




IRAM 


1-0 


<2.6 


0.3 


2 


UCyg 


OSO 


1-0 


4.9 


4.8 


2 




IRAM 


2-1 


2.6 


2.8 


2 




IRAM 


2-1 


9.4 


31.7 


2 














1. This paper; 2. Olofsson et al. ( 


1993a 


); 3. Neri et al. ( 


1998); 4. JCMT public archive. 



1.5 x 10 _5 Mq yr -1 , over which the physical conditions of 
the CSEs vary considerably, posing a challenge to the model. 
The intensities, overall line shapes, and when available radial 
brightness distributions, of the circumstellar CO lines produced 
by the radiative transfer model generally agree well with those 
observed. This is illustrated here for three of our sample stars 
LP And (Fig. @), R For (Fig. g), and W Ori (Fig. §). These 
carbon stars span a large range in mass loss rate and serve to 
illustrate the various physical conditions in these CSEs. 

A full, detailed, error analysis of the estimated mass loss 
rates is not possible due to the relatively large number of free 
parameters entering the model. Instead, we will vary some of 



the more important parameters in order to illustrate the sensi- 
tivity of the model, and to be able to get a rough estimate of the 
errors involved in the mass loss rate estimates. The results of 
these sensitivity tests are shown in Table gfor our three exam- 
ple stars. In addition, in Fig. ^ we present chi-square contour 
plots for these stars, produced by varying the mass loss rate 
and the /i-parameter. This will give an estimate of the accu- 
racy in the determination of these two adjustable parameters 
when all other parameters are held fixed (the size and shape of 
the CO envelope is however allowed to change in accordance 
with Eqs [l(] and |TT]). In Fig. M we indicate various confidence 
levels with the contour at x^~ n +2.3 marking the 68% confi- 
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Fig. 4. Multi-transition CO millimetre-wave line emission observed towards the high mass loss rate Mira variable LP And. 
The observed spectra (histograms) have been overlayed with the model prediction (full line) using a mass loss rate of 
1.5x 10 _5 Mq yr _1 . The transition, telescope used, and the corresponding beamsize, are indicated for each of the observations. 
Cooling panel: The full line represents adiabatic cooling; the dotted line gives the H2 cooling; and the dashed line is CO cool- 
ing (see text for details). Temperature/optical-depth panel: The dotted line shows the kinetic gas temperature derived from the 
energy balance equation. The full line gives the excitation temperature of the CO( J=2^1) transition. The dashed lines give the 
tangential optical depth, T^ an , of the CO( J=2^1) transition. 



dence limit (i.e., the "la" level for two adjustable parameters). 
The information contained in the shape of the line profiles (to 
which the \ 2 defined above is not very sensitive) have occa- 
sionally been used to further constrain the most probable pa- 
rameter space. 

LP And is a high mass loss rate Mira variable where the ex- 
citation of 12 CO is dominated by collisions. Consequently, the 



line intensities are very sensitive to the temperature structure, 
i.e., the ^.-parameter. It is interesting to note that the assumed 
envelope size r p does not significanly affect the derived mass 
loss rate in the high mass loss rate regime. The reason for this 
being that the density is too low to excite the CO molecules ef- 
fectively (the lines are sub-thermally excited) in the cool outer 
parts of the envelope (cf. Fig. 0), i.e., the emission is excitation 
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Fig. 5. Multi-transition CO millimetre-wave line emission observed towards the bright carbon star R For. The observed spec- 
tra (histograms) have been overlayed with the model prediction (full line) using a mass loss rate of 1.3x 10~ 6 M Q yr _1 . The 
transition, telescope used, and the corresponding beamsize, are indicated for each of the observations. Cooling panel: The full 
line represents adiabatic cooling; the dotted line gives the H2 cooling; and the dashed line is CO cooling (see text for details). 
Temperature/optical-depth panel: The dotted line shows the kinetic gas temperature derived from the energy balance equation. 
The full line gives the excitation temperature of the CO(J=2^1) transition, and a dash-dot line indicates a negative excitation 
temperature, i.e., maser action. The dashed lines give the tangential optical depth, r tan , of the CO( J=2— >1) transition. 



limited. In the high mass loss rate regime the line intensities 
are also insensitive to the adopted mass loss rate since an in- 
crease in M leads to more cooling (CocM while iJcxM 5 ), 
which compensates for the increase of molecular density. This 
"saturation"-effect of the line intensities for high mass loss 



rates has been noted before in other models where the cool- 
ing by CO is treated in a self-consistent manner (Sahai 1990; 
Kastner 1992). Additional sensitivity tests in the case of a high 
mass loss rate object (CW Leo) can be found in Groenewegen 
(1994). 
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Table 5. The effect on the integrated model intensity 7 mo d (in percent), obtained by SEST and JCMT, due to changes in various 
parameters. In general, the accuracy of th observed intensities lie at the ±20% level. 



LP And R For W Ori 



Parameter 


Change 


1-0 


2-1 


3-2 


4-3 


1-0 


2-1 


3-2 


4-3 


1-0 


2-1 


3-2 


4-3 


M 1 ' 2 


+50% 


+25 


+15 


+ 15 


+10 


+50 


+55 


+55 


+55 


-10 


+ 15 


+40 


+55 




-33% 


-25 


-15 


-15 


-15 


-35 


-35 


-40 


-40 


-40 


-30 


-25 


-35 


M 


+50% 


+ 15 








+5 


+55 


+40 


+35 


+30 


+25 


+30 


+40 


+45 




-33% 


-15 


-5 


-5 


-5 


-35 


-30 


-30 


-25 


-40 


-25 


-20 


-30 


L 1 


+50% 


























+ 15 


+ 10 


+5 


+5 




-33% 

















-5 








-20 


-10 


-5 


-5 


h 


+50% 


+25 


+35 


+35 


+40 


+ 10 


+ 15 


+20 


+25 





+5 


+ 10 


+ 15 




-33% 


-20 


-25 


-30 


-30 


-10 


-15 


-20 


-20 





-10 


-10 


-15 


r P 


+50% 


+5 











+ 15 


+5 


+5 





+20 


+5 


+5 







-33% 


-5 











-20 


-5 








-50 


-30 


-10 


-5 


n 


x2 


+5 











+5 


+5 


+5 





+5 


+ 15 


+ 15 


+ 10 


Vt 


x2 


+ 10 











+ 15 


+ 10 


+5 


+5 





+25 


+40 


+30 



1 Same temperature structure as final model was used; 2 Same envelope size as final model was used. 



W Ori, on the other hand, is a low mass loss rate object 
with low optical depths in the lowest rotational transitions. In 
fact, some of the low lying transitions are inverted over parts 
of the CSE, i.e., maser action is taking place. In particular, the 
J=1^0 is inverted over the entire CSE producing a line pro- 
file that is clearly double peaked, i.e., this is due to maser action 
and not due to resolution effects. The observations of this tran- 
sition are complicated by possible contamination by interstel- 
lar lines, and the signal-to-noise ratio is not particularly good, 
erasing any possible signs of maser action. As expected the ra- 
diation emitted by the central star plays an important role in the 
excitation. The transitions exhibiting population inversion are 
particularly sensitive (both the intensity and line profile) to the 
physical conditions prevailing in the CSE. In this case where 
radiative excitation is important the model becomes sensitive 
to the choice of inner radius, r- u as well as tubulent velocity, v t . 
Here, also the size of the envelope is important when determin- 
ing the mass loss rate, i.e., the emission is photodissociation 
limited (at least for the lowest transitions). Generally, the line 
intensities scale roughly linearly with M. 

R For presents an intermediate case between LP And and 
W Ori. It is also a typical sample star in the sense that the ma- 
jority of the stars have envelope parameters that resemble those 
of R For. From Table || and Fig. [7] we see that R For shares prop- 
erties of both the low and the high mass loss rate star. An inter- 
mediate mass loss rate model is generally more sensitive to the 
temperature structure than the radiation field, thus resembling 
the high mass loss rate objects. However, as in the low mass 
loss rate regime the line intensities scale roughly linearly with 
M. We also note that the derived line intensities, especially the 
J=1^0 transition, start to become sensitive to the outer ra- 
dius of the CO envelope. The exact mass loss rate where the 
transition from photodissociation- to excitation-limited emis- 
sion occurs is hard to pinpoint, since it depends on the envelope 
characteristics. We have made tests which show that for typical 



envelope parameters the transition lies at ^5xlO _7 M Q yr _1 
for the J=1^0 transition. Higher J-lines sample hotter and 
denser gas, located closer to the star, and are therefore less sen- 
sitive to the choice of the envelope size. 

In summary, we believe that within the adopted circumstel- 
lar model, the estimated mass loss rates are accurate to about 
±50% (neglecting errors introduced by the uncertain CO abun- 
dance and the distance estimates) when good observational 
contraints are available, but one should keep in mind that the 
causes of the uncertainty varies with the mass loss rate. 

Another rough way to estimate the errors involved is to 
compare with the results obtained from other self-consistent 
models. When comparing our mass loss rate estimate for CW 
Leo to those obtained from detailed radiative transfer models 
(Kastner 1992; Crosas & Menten 1997; Groenewegen et al. 
1998a; Skinner et al. 1999) we find a very good agreement, 
within 20 %, when adjustments for differences in fo and dis- 
tance have been made. Kastner (1992) also modelled the high 
mass loss rate object RW LMi (a.k.a. CIT 6) and obtained a 
mass loss rate in excellent agreement with our estimate (when 
corrected for the difference in distance). Recently, Hiriart & 



Kwan (2000) presented a method to combine a circumstellar 



dust model with a CO model, keeping physical quantities con- 
sistent between the models. They present mass loss rate esti- 
mates for 17 of our sample stars. Based upon their Eq. 17 and 
Table 6 we derive hydrogen mass loss rates for these stars that 
are, on the average, 1.3 times larger than those derived from 
our radiative transfer analysis (for an individual star, however, 
the discrepancy can be as large as a factor of 3). 

We also compare our derived mass loss rates with those 
obtained by Olofsson et al. (1 993a) using a formula based on 
the work by Knapp & Morris (1985), 



Mi 



KM 



= 5.7 x 10 



Tn,uV„D 



-20 mb 



2fl2 



nib 



s(J)fS 



0.85 



M yr" 



(14) 
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Fig. 8. The ratio of the estimated mass loss rates from our radia- 
tive transfer analysis (M) to those estimated from a analytical 
formula (Mrm) (see text for details). 



which relates the mass loss rate Mkm to easily determined ob- 
servables (T m b is given in Kelvins, v c in kms -1 , D in pc, and 
6* m b in arcseconds). In Eq. [14] s( J) is a factor that depends on 
the transition (J) used. For J=1^0 s(J)=l and for J=2^1 
s( J)=0.5 (see the discussion in Olofsson et al. 1993a). When 
estimating the mass loss rates from Eq. [l4| we have used the ob- 
served J=l— >0 and J=2^1 line emission, and averaged the 
results. This mass loss rate was compared to that obtained from 
the radiative transfer analysis, Fig. |[ It is found that Eq. |l4| sig- 
nificantly underestimates the mass loss rates when compared 
to those obtained from the radiative transfer analysis by, on 
the average, a factor of about four. However, considering the 
simplicity of Eq. |l4| the agreement is remarkably good. There 
is also a trend that the discrepancy increases with lower mass 
loss rate (indeed Olofsson et al. suspected that their simple ap- 
proach led to a systematic underestimate of the mass loss rate 
that became worse the lower the mass loss rate). In this context 



it should be noted that Eq. 14 was derived assuming a fixed 
CO envelope size of 3 x 10 17 cm, which is appropriate for high 
mass loss rate objects. In the high mass loss rate regime the 
line intensities are not very sensitive to the exact envelope size 
since the emission is excitation limited. This is, however, gen- 
erally not the case for the optically bright carbon stars in our 
sample, where the line intensities are sensitive to the size of the 
CO envelope set by photodissociation. A correction for the en- 
velope size [see Neri et al (1998) and references therein] would 
decrease the discrepancy seen in Fig. ^|for the stars losing mass 
at lower rates. Taken at face values the discrepancies between 
the mass loss rates derived here and those derived by Olofsson 
et al. (1993a) are smaller, but this is due to the fact that on the 
average Olofsson et al. used larger distances. This also means 



that the trend in the dust-to-gas mass ratio, decreasing with the 
mass loss rate, that Olofsson et al. reported was due to underes- 
timated gas mass loss rates [see also Hiriart & Kwan (2000)]. 



Indeed, in Sec. 1.2 we find the opposite behaviour, i.e., indica- 



tions of a dust-to-gas mass ratio that increases with the mass 
loss rate. 



4.4. The CO envelope sizes 

The size of the CO envelope is an important parameter, and 
the derived mass loss rate will depend on it to an extent that 
depends on the mass loss rate (see discussion in Sec. 4.3). 
Neri et al. ( |1998| ) presented IRAM PdB interferometer (PDBI) 
CO(J=1^0) maps of varying quality for 9 of our sample 
stars. From their CLEANed brightness maps we have obtained 
the averaged radial brightness distributions. In addition, we 
have obtained two single-dish CO(J=1^0) maps using the 
OSO 20 m telescope. In Fig. ^ we compare the observed radial 
brightness distributions to those calculated from our model, 
which uses the Mamon et al 



results to estimate the 
CO envelope size, for the parameters presented in Tab. [3| We 
find that the estimated CO envelope sizes are generally consis- 
tent with the observations. In four cases (V384 Per, RW LMi, 
T Dra, and S Cep) the calculated radial brightness distribu- 
tions fail to reproduce the observations. At least in the case 
of V384 Per and RW LMi this appears to be due to deviations 
from a simple r~ 2 density law, in the sense that two distinct 
components are clearly visible in the brightness maps. 

We are not able to distinguish between the different CO 



envelope sizes obtained by Mamon et al. (1988) and Doty & 



Leung (1998 ) for the high mass loss rate object CW Leo. This 



is simply because the emission in this case is excitation limited 
with the external cool and tenuous parts producing very little 



emission. 



4.5. The kinetic temperature structure 



As was demonstrated in Sec. 1.3 the derived mass loss rate can 



be very sensitive to the kinetic temperature structure of the gas. 
Adiabatic cooling alone results in T(r)ocr s , where S=— 4/3 
(assuming 7=5/3), but from Figs. ^|it is evident that other 
cooling processes as well as heating mechanisms are important. 
If CO line radiation is the dominant cooling mechanism and 
n>n CI i t one would expect T(r)ocr _1 with a somewhat shal- 
lower decline if n<n CI ;t (see Doty & Leung 1997). 

Generally, it is not possible to assign a simple power law 
to the radial behaviour of the temperature throughout the en- 
tire CSE. However, in the case of LP And we find 5^—0.9 for 
r<5xl0 16 cm. Further out in the cool, tenuous, outer layers 
of this CSE, heating due to the photoelectric effect is impor- 
tant and raises the temperature. For R For and W Ori we esti- 
mate, over the regions where the observed emissions emanate, 
i.e., from — lxlO 15 to — 5xl0 16 cm, S to be —0.7 for both 
sources. For comparison, Doty & Leung (1997 ) derive 1 .0 



for the high mass loss rate object CW Leo, while we obtain 
S — 0.9. 
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From the fact that CocM and HocM 5 one would expect 
high mass loss rate objects to have significantly cooler en- 
velopes than low mass loss rate objects (e.g., Jura et al. 1988). 
In our modelling we do not see this trend. The reason being that 
in addition Hoth(Lv c ) 3 / 2 , where all of the parameters have a 
mass loss rate dependence in the sense that their values are gen- 
erally smaller for the low mass loss rate objects. 

4.6. 13 CO line cooling 

Carbon stars generally have 12 C/ 13 C-ratios in the range 30—70 
(Lambert et al. 1986; Schoier & Olofsson 2000). This means 
that for the large majority of all stars in our sample line cooling 
by 13 CO is insignificant, even when optical depth effects are 
taken into account. However, there exists a small population of 
carbon stars, the J-type stars, which have 12 C/ 13 C-ratios of ~3. 
For these stars 13 CO line cooling could be important. 

To test this we have included 13 CO line cooling for one of 
our sample stars, Y CVn, a known J-type star. We h ave u sed 
the 13 CO model as presented in Schoier & Olofsson (2000) to 
calculate the the contribution from this isotopomer to the total 
amount of cooling, and used this extra cooling term as input 
to the standard model. We find that, after we have iterated be- 
tween the 12 CO and 13 CO models a few times, that 13 CO line 
cooling is somewhat less effective as a coolant than expected, 
based on the isotope ratio (probably due to slightly different 
excitation conditions between the two molecules). The effect 
on the derived line intensities is small; the largest effect is on 
the J=3— >2 line emission which is lowered by ^10%. The 
changes in line intensity may be compensated for by increasing 
the ^.-parameter. 



Moreover, Ryde et al. ( 1999 ) present a model of the high 
mass loss rate object IRAS+15 194-51 15 using the radiative 
transfer code outlined here. The 12 C/ 13 C-ratio for this object 
was estimated to be 5.5, i.e., it resembles the J-type stars in this 
respect. It was found that 13 CO line cooling in this case was 
approximately one third of that of 12 CO (due to optical depth 
effects the 13 CO line cooling is in a relative sense more impor- 
tant). The decrease in kinetic temperature due to this extra cool- 
ing was compensated for by the fact that the /i-parameter (heat- 
ing) had to be somewhat increased to produce the observed line 
intensities. Consequently, the derived mass loss rate was not af- 
fected by introducing the 13 CO line cooling. 

5. Deviations from the standard model 

For some stars, where there are enough observational con- 
straints, it is possible to establish deviations from the standard 
model. 



5.7. Detached CSEs 

In our sample there exist a number of stars that exhibit a 60 /im 
excess, placing them in the regions Via and VIb in the IRAS 
colour-colour diagram (see Olofsson et al. 1993a), suggesting 
that their CSEs are abnormal in some way. Previous CO obser- 



vations of five of these 60 /im excess stars (see Table ^) have 
revealed that their CSEs cannot have been produced by a single 
smooth expanding wind, but rather suggested a drastic modu- 
lation of the mass loss on a time scale as short as a thousand 
years, which produced detached CSEs (dCSEs; Olofsson et al. 
1996b). Presently, these stars are losing matter at a significantly 
lower rate, forming what is refered to as an attached circumstel- 
lar envelope (aCSE). 

To be able to put constraints on the mass loss and its varia- 
tion with time, high spatial resolution observations are needed. 
Such observations have been carried out using the IRAM PdB 
interferometer, with an angular resolution of for two of 
the stars in Table |(| U Cam and TT Cyg. It was found that the 
dCSE around U Cam is fairly young while that around TT Cyg 
is significantly older. Both dCSEs show a remarkable overall 
spherical symmetry, and they are geometrically very thin. The 
CO line emission towards these two sources, supplemented by 
a detailed modelling using the radiative transfer code outlined 



in this paper, were presented in Lindqvist et al. (1999) (U Cam) 
and Olofsson et al. (1998, 2000) (TT Cyg). 

In this paper we have attempted to model the remaining 
three stars with known detached shells (R Scl, S Set, and 
U Ant), based upon the results obtained from the high angular 
resolution observations described above. The results for both 
the aCSEs and dCSEs around these stars are presented in Ta- 
ble H (where also the results for TT Cyg and U Cam are in- 
cluded). Due to the limited observational constraints available, 
the dCSEs around S Set and U Ant were modelled using the 
same geometrical thickness and kinetic temperature structure 
as obtained for TT Cyg (Olofsson et al. 2000). For R Scl we 
used the same geometrical thickness as adopted by Olofsson 
et al. 1996b and assumed a kinetic temperature decrasing lin- 
early with radius from a value of 50 K at the inner shell radius. 
The derived intensities from the model generally agree well 
with those observed (Table ^|). However, we note that in the 
case of R Scl we only get about half of the observed J=3— >2 
line emission from the model, but here the aCSE could pro- 
duce a significant amount of the observed emission. We find 
that the estimated shell masses, M s , are relatively similar, but 
there is a trend of increasing shell mass with age (the R Scl 
shell mass may be overestimated due to the difficulty in sepa- 
rating the aCSE and dCSE emission), which may indicate that 
gas is being swept up. 

The present mass loss rate was modelled for U Ant and 
S Set using the standard model. The results are presented in 
Table ^ and Table [7| For R Scl the emission from the aCSE 
and dCSE cannot be separated. Again, the predicted line in- 
tensities agree, with one exception, well with those observed, 
considering the uncertainties of the observations. For S Set the 
model of the aCSE fails to reproduce the IRAM CO( J=2— >1) 
observations, in contrast to the dCSE model where the fit to the 
observations is excellent. An explanation to this inconsistency 
could be that there are pointing problems in these observations 
(the IRAM observations are more sensitive to the pointing due 
to its smaller beam), which would mainly affect the emission 
observed from the aCSE. We find that the present mass loss 
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aCSE 



dCSE 



Source 


Var. type 


P 


D 


M 






M B 






r a 


Ar 8 


M s 








[days] 


[pc] 


[Moyr- 1 ] 


[kms^ 1 ] 


[cm] 


[M yr" 


'] 


[kms -1 ] 


[cm] 


[cm] 


[M ] 




RScl 1 


SRb 


370 


360 








4.0x10" 


-6 


16.5 


7.0xl0 16 


7.0xl0 16 


5.4x10" 


-3 


UCam 2 


SRb 




500 


2.5xl0~ 7 


12.0 


1.4xl0 16 


9.0x10" 


-6 


23.0 


5.5xl0 16 


l.lxlO 16 


1.4x10" 


-3 


U Ant 1 


Lb 




260 


3.0xl0~ 8 


4.5 


1.4xl0 16 


1.5x10" 


-5 


19.5 


1.6xl0 17 


2.0xl0 16 


4.9x10" 


-3 


S Set 1 


SRb 


148 


400 


3.0xl0- 8 


4.5 


1.4xl0 16 


4.0x10" 


-5 


17.0 


4.0xl0 17 


2.0xl0 16 


1.5x10" 


-2 


TTCyg 3 


SRb 


118 


510 


3.0xl0~ 8 


4.0 


1.5xl0 16 


1.5x10" 


-5 


12.5 


2.7xl0 17 


1.9xl0 16 


7.2x10" 


-3 



Model is presented in: 1 This paper; 2 Lindqvist et al. (19 9p); " p iolsson et al. (2000 | | 

Table 7. CO modelling results, for three sources with known detached shells, compared to observations. 



aCSE 



dCSE 



Source 


Tel. 


Trans. 


-^obs 


-^mo d 


-^obs 


-^mod 


Rel. 








[Kkms -1 ] 


[Kkms -1 ] 


[Kkms -1 ] 


[Kkms -1 ] 




RScf 


SEST 


1-0 






26.1 


29.6 


2 




IRAM 


1-0 






76.2 


79.2 


1 




SEST 


2-1 






47.8 


49.8 


1 




JCMT 


2-1 






50.6 


54.8 


3 




IRAM 


2-1 






67.7 


72.0 


1 




SEST 


3-2 






63.7 


28.0 


2 




JCMT 


3-2 






62.9 


30.5 


3 


U Ant 


SEST 


1-0 




0.6 


10.0 


13.4 


2 




SEST 


2-1 


3.2 


2.8 


7.6 


7.3 


2 




JCMT 


2-1 


3.6 


3.3 


7.7 


7.2 


3 




SEST 


3-2 


3.6 


4.0 


5.0 


3.0 


2 


S Set 


SEST 


1-0 




0.3 


5.3 


5.6 


1 




IRAM 


1-0 


2.5: 


1.3 


5.3 


5.4 


1 




SEST 


2-1 


1.1 


1.3 


2.7 


2.9 


2 




JCMT 


2-1 


0.9 


1.3 


2.4 


3.0 


3 




IRAM 


2-1 


1.8: 


5.4 


2.8 


2.8 


1 




SEST 


3-2 


2.1 


1.7 


1.0 


0.9 


2 




JCMT 


3-2 


2.4 


2.2 


1.0 


0.9 


3 



a For this source it is not possible to separate the contributions from an aCSE and a dCSE; 1. Olofsson et al. ( 1993a); 2. Olofsson et al. ( 1996) 
3. JCMT public archive 



rates of the stars with the oldest shells are low, and that the 
present mass loss rate of U Cam (which has a much younger 
dCSE) is significantly higher. 

The central stars are semi-regular (SR) or irregular vari- 
ables (Lb), and they have low present mass loss rates. They ap- 
pear to have been subject to drastic modulations of their mass 
loss rates over relatively short periods of time (a few thousand 
years). We estimate that the CO dCSE emission is observable 
for ~10 4 years (Bergman et al. 1993). If this mass ejection is 
a repeatable phenomenon, and all our sample stars go through 
this phase, we estimate that the time scale between mass ejec- 
tions is about 10 5 years. Therefore, the mass ejections may pos- 
sibly be linked with the helium-shell flashes (the process that 
leads to nuclear-processed matter being transported to the sur- 
face of the star) predicted to occur regularly in AGB stars with 
time intervals of the order of 10 5 years for the relevant masses 
(e.g., Blocker 1995). We note that all carbon stars that have 
been found to have clearly detached gas shells are members of 



our sample, and that no M-type star with a detached CO shell 
has been detected. 



5.2. V Hya and TX Psc 

The sucessful modelling of the vast majority of our sample 
stars, assuming a spherically symmetric CSE, suggests that 
axi- or non-symmetric mass loss is not a common phenomenon 
among these optically bright carbon stars. However, in the case 
of V Hya and TX Psc, stars with complex radio line profiles, 
the possibility of bipolar outflows has been considered (Heske 
et al. 1989; Kahane et al. 1996). Due to the complexity of these 
outflows, no radiative transfer analysis was attempted here. 



Kahane et al. (1996) suggest that V Hya, in addition to the 
high velocity bipolar outflow, has a slowly expanding spheri- 
cal component. Using a 2D radiative transfer code, based upon 
the Sobolev approximation, they derive a total mass loss rate of 
~1.5x 10 _6 M Q yr _1 for this object. The authors further sug- 
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gest that V Hya is a transition object, between carbon stars and 
PNe, having just developed a highly axi-symmetric mass loss. 
The slowly expanding spherical component is then interpreted 
as the fossil CSE created during the AGB stage. Knapp et al. 



10" 



(1999), however, argue, based upon the spectral type, period, 



colours, and lack of ionizing radiation, that V Hya is still on 
the AGB. Moreover, they estimate the total mass loss rate of 
this object to be ^4x 10 _5 M Q yr~ 1 . This, in addition to the 
small dynamical age of the envelope, suggests that V Hya has 
entered its "superwind" phase. 

The observed line profiles of the CO line emission towards 
TX Psc are somewhat peculiar (Heske et al. 1989; Olofsson 
et al. 1993a). As a consequence, the expansion velocity de- 
rived from the observed radio line profiles range from 7.5 to 
12.2 km s -1 (see Tableland Olofsson et al. 1993a) depending 
on the transition and the telescope used. Heske et al. (1989), 
mapped the circumstellar CO line emission around this object, 
and interpreted the observed asymmetries as being produced ei- 
ther by a bipolar outflow or a highly clumped wind. No detailed 
modelling exist for the CSE around TX Psc, and its evolution- 
ary status remains uncertain. It might be in a state similar to 
that of V Hya. 

5.3. DR Ser 

DR Ser has the highest 60 /im excess of all the sample stars, 
which places it in the IRAS two-colour diagram among objects 
having detached envelopes (see Table |6|). Indeed, Kerschbaum 
( 1999 ) suggests that this might be a detached shell source based 



upon his modelling of the SED of this object. However, the CO 
single dish observations reveal no signs of double-peaked pro- 
files (confusion with interstellar lines complicates the interpre- 
tation), which are indicative of a detached CO shell. Moreover, 
DR Ser has a relatively low 12 C/ 13 C-ratio of 6 (Abia & Isern 
1997). Thus, it resembles the J-type stars. 

Attempting to model this source with a standard model of 
its CSE fails to explain both the observed J=1^0 and J=2^1 
line intensities in that they come out too weak, by more than a 
factor of two. This is due to the extremely cold envelope, result- 
ing from the adopted h -parameter of 0.2 (due to the luminosity 
being lower than 6000 Lq), in combination with a relatively 
high mass loss rate ~2x 10~ 6 M Q yr~ 1 required to maximize 
the line intensities. Further increasing the mass loss rate will 
make the envelope even cooler due to the increase in CO line 
cooling. If one instead assigns an /i-parameter of ^1 it is pos- 
sible to obtain a reasonable fit to the data using a mass loss 
rate of ^5x 10 -6 M Q yr -1 . This estimate must be regarded as 
highly uncertain. 



6. Discussion 

6.1. The mass loss rate distribution 

Fig. [T(^ shows the distribution of the mass loss rates obtained 
from the radiative transfer model for all sources (also shown is 
the distribution for the stars within 500 pc of the Sun), as well 
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Fig. 11. The derived mass loss rate plotted against the gas ex- 
pansion velocity of the CSE. Sources with known detached 
shells are indicated with triangles (open triangles represent dC- 
SEs, while filled ones represent aCSEs). M scales as wf , with 
£=2.5. Also shown is the correlation coefficient r. 



as divided into variability groups (note that objects with known 
detached shells have not been included in this analysis). Since 
the statistics of the SRa stars are poor we have chosen to group 
them together with the SRb stars. The mass loss rate distribu- 
tion for all stars is sharply peaked around the median mass loss 
rate of 2.8 x 10 _7 Mq yr _1 (the median mass loss rate for the 
stars within 500 pc of the Sun is 1.6x 10 _7 Mq yr _1 ). It is in- 
teresting to note that this is roughly the rate at which the core 
mass is expected to grow due to nuclear burning (Schonberner 
1983). Mira variables generally have larger mass loss rates than 
other variability types, while irregular variables and semiregu- 
lars appear to have very similar mass loss rate characteristics. 

Based on the mass loss rate distribution for the sample 
of stars within 500 pc, which we believe to be close to com- 
plete and for which we have detected all sources, we are able 
to draw some general conclusions. The sharp decline at mass 
loss rates below ^5x 10 _8 M Q yr~ 1 is very likely real. Netzer 
& Elitzur (1993j) estimate that a mass loss rate in excess of 



~10 -7 Mq yr _i is required to get a dust-driven wind. The ex- 
act limit is, however, sensitive to the adopted stellar and dust 
parameters. The drastically decreasing number of high mass 
loss rate objects is also real, although our selection criterion 
bias against these objects, and can be explained by the fact that 
carbon stars of the type discussed here, i.e., mainly low mass 
ones (Claussen et al. 1987), only for a limited time, or possibly 
never, reach high mass loss rates. 
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6.2. Mass loss and envelope kinematics 

The v c -distribution for this sample of stars has already been 
shown and discussed by Olofsson et al. ( 1993a| ). However, 
a new comparison between the mass loss characteristics M 
and v e is warranted considering the more reliable mass loss 
rates obtained in this paper, Fig. O. We find a clear trend 
that v c increases with M, v e <xM° (with a correlation co- 



efficient of 0.78). In comparison, Olofsson et al. (1993c) de- 
rived v e (xM ' 53 , but their mass loss rates were calibrated us- 
ing Eq. [l4], which includes a v% -dependence. The correlation 
is much tighter than that obtained by Olofsson el al. (1993a), 
which is reassuring. Thus, the mass loss mechanism operates 
such that mass loss rate and expansion velocity increase to- 
gether. However, the scatter appears larger than the uncer- 
tainties in the estimates, and so the mechanism also produces 
widely different mass loss rates for a given expansion veloc- 
ity. At high mass loss rates the scatter is larger with a possible 
divison into objects with high mass loss rates but only mod- 
erately high velocities, and objects with moderately high mass 
loss rates but high velocities. We note also that for the detached 
shell sources the mass loss rates and expansion velocities for 
the two mass loss epochs follow the general trend. 



Habing et al. (1994) studied the momentum transfer from 
the photons, via the dust, to the gas in the CSEs of AGB stars. 
They found that for mass loss rates close to the minimum 
mass loss rate for a dust-driven wind, ^3x 10 -8 Mq yr _1 , v c 
increases linearly with M, while for mass loss rates above 
~10 -6 Mq yr _1 the dependence weakens considerably. At 
10~ 5 MQyr -1 they found v e ocM 0M . Habing et al. attribute 
the increase in v c with mass loss rate (for low mass loss rates) 
to a higher efficiency in the coupling between gas and dust. 
This can explain the observed behaviour in Fig. [n] for mass 
loss rates below about 10 -6 Mq yr _1 , but for the higher mass 
loss rates it appears that a dependence of v c (and M) on lumi- 



nosity gives the most reasonable explanation, see Sect. 6.3 



6.3. Dependence on stellar properties 

In Fig. [12] we plot the circumstellar characteristics, M and v c , 
against the stellar characteristics luminosity (L), period (P), ef- 
fective temperature (T ff ), and the photospheric C/O-ratio. The 
effective temperatures and C/O-ra tios use d in this analysis are 
those presented in Olofsson et al. ( 1993b ). We have looked for 
dependences assuming that the ordinate scales as the abscissa 
to the power of £, and a normal correlation coefficient r was 
calculated as an estimate of the quality of the fit, see Fig. [l~2|. 
The uncertainties in the estimated quantities are of the order: 
±50% (M), ±2kms- 1 (v e ), ±10% (P; but some periods may 
be poorly determined), a factor of 2 (L), and ±200 K (T e g). 

Clearly, both the mass loss rate and the expansion velocity 
increase with the pulsation period of the star. There is a weaker 
trend with the luminosity of the star. The latter dependence is 
not completely independent of the former, since some of the lu- 
minosities are estimated from a P — L relation. The existence 
of a P — L relation is usually attributed to a distribution in mass 



(Jones et al. 1994), i.e., the higher the mass the longer the pe- 
riod and the higher the luminosity. When looked at in detail, at 
least the apparent mass loss rate dependence on period may be 
attributed to a change from semiregular pulsation at short pe- 
riods to regular pulsations at longer periods. We may therefore 
infer that for the mass lass rate it is not clear whether it is the 
regularity of the pulsation or the luminosity that causes the in- 
crease with period. For the expansion velocity the increase with 
period may be a combined effect of increasing mass loss rate 
with period (for low mass loss rates) and an incr ease i n lumi- 
nosity for the longer periods, e.g., Habing et al. ( 1994 ) derive 
v e ocL - 35 in their dust-driven wind model. In addition, there is 
evidence of a weak trend in the sense that higher mass loss rate 
objects have lower effective temperatures. These correlations 
are all consistent with a dust-driven wind, where the pulsation 
may play an important role. There appears to be no correlation 
between the mass loss rate and the C/O-ratio, which is surpris- 
ing considering that the dust-to-gas mass ratio in a C-rich CSE 
should be sensitively dependent on this. 

The present mass loss rates for the detached shell sources 
are typical for their periods but somewhat low for their lumi- 
nosities. However, during the formation of the dCSEs the mass 
loss rate must have been atypically high for the present peri- 
ods and luminosities. We also find that the expansion velocities 
of the dCSEs lie at the very high end of expansion velocities 
found for other stars with the same luminosity, indicating that 
these stars had higher luminosities during the shell ejection. 

Two other stars, SZ Car and WZ Cas, stand out in these 
plots. The properties of SZ Car resembles the stars with known 
dCSEs, but the CO line profiles give no indication of a detached 
shell. WZ Cas has by far the lowest estimated mass loss rate in 
the sample, as well as the lowest C/O-ratio, 1.01, which clas- 
sifies it as an SC-star. It is also a Li-rich J-star (Abia & Isern 
1997). 

6.4. Enrichment of the ISM 

Carbon stars on the AGB are important in returning processed 
gas to the interstellar medium (ISM). The total mass loss rate 
of carbon stars in the Galaxy is obtained from 



M Gal = 



2irRi:<M>dR 



(15) 



where £ is the surface density of carbon stars, <M> is 
the mean mass loss rate of the carbon stars. According to 
Guglielmo et al. (1998) the infrared carbon stars have a roughly 
constant surface density out to the galactocentric distance of the 
Sun, after which it follows an exponential decline with a scale 
length of ^2.5 kpc. We approximate this by assuming a con- 
stant value, So (the surface density in the solar neighborhood), 
to lOkpc and zero beyond this galactocentric distance. 

Based on the complete sample, i.e., stars within 
500 pc from the Sun, we estimate Y>o<M> to be 
^1.7x 10~ 10 Mq yr _1 pc~ 2 , where we have included 
also helium. This estimate is very sensitive to the number of 
high mass loss rate objects found within 500 pc, e.g., the high 
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mass loss rate object CW Leo contributes almost half of our 
estimate of the total mass returned to the ISM. Our estimate 
of the rate at which matter is returned to the ISM by carbon 
stars is consistent with previous estimates (considering the 
large uncertainties), e.g., Knapp & Morris (1985) and Jura & 

Mq yr -1 pc~ 2 



-To" 



Kleinmann (1989|) derive values of ~2x 10 

2 respectively. Using Eq. |15| 



and 



4.5x10^™ M yr _1 pc 



the annual return of matter to the ISM in the Galaxy is 
estimated to be ^0.05 Mq for the carbon stars considered 
here. It is quite possible that a larger mass return from carbon 
stars is obtained during their final evolution on the AGB, a so 
called superwind phase, as is indicated by the estimated values 
for infrared carbon stars, ~0.5 Mq yr" 1 (Epchtein et al. 1990; 
Guglielmo et al. 1993 ), and PNe, ~O.3M yr" 1 (Gustafsson 
et al. 1999). This confirms the importance of carbon stars 
for the cosmic gas cycle in galaxies. By comparison, high 
mass stars are estimated to contribute about 0.04MQyr _1 
(Gustafsson et al. 1999). 

We estimate the contribution from carbon stars to the car- 
bon enrichment of the ISM to be ~0.5x 10 _4 Mq yr" 1 during 
the AGB stage (if the subsequent superwind phase is included 
this value may increase to ~5x 10 _4 Mq yr" 1 ). This corrobo- 



rates the conclusions by Gustafsson et al. ( 1999 ) that 'normal' 
carbon stars are not important for the total carbon budget of 
the Galaxy. According to these authors the main contributors 
should instead be high mass stars in the Wolf-Rayet stage, an- 
nually supplying the Galaxy with ^0.01 Mq of carbon. This is 
roughly what is required to produce the ^10 _3 x10 11 Mq of 
carbon present in the Milky Way over a period of 10 10 years. 



7. Conclusions 

We have developed a detailed radiative transfer code to model 
circumstellar molecular line emission. The code also solves for 
the energy balance equation of the gas. It is found that the mass 
loss rate determination for low mass loss rate objects depends 
crucially on a number of assumptions in the CSE model, ex- 
cept on the dust properties, since the CO molecules are radia- 
tively excited. On the other hand, the mass loss rate estimate 
for a high mass loss rate object depends essentially only on the 
temperature structure, and hence the uncertain dust parameters. 
In addition, for such an object the CO emission saturates, and 
becomes less useful as a mass loss rate measure. We also find 
that different lines respond differently to changes in the various 
parameters. Therefore, a reliable mass loss rate determination 
requires, in addition to a detailed radiative transfer analysis, 
good observational constraints in the form of multi-transition 
observations and radial brightness distributions. 

This model has been applied to CO radio line observations 
of a large sample of optically bright N-type and J-type carbon 
stars on the AGB (69 objects). The sample is reasonably com- 
plete out to about 1 kpc, and all stars (41) within ^500 pc of 
the Sun have been detected in circumstellar CO line emission. 
The derived mass loss rates span almost four orders of mag- 
nitude, — 5xlO" 9 M Q yr" 1 to — 2x 10" 5 M Q yr _1 , over which 
the physical conditions of the CSEs vary considerably. The fact 



that the model can be succesfully applied over such a wide 
range of environments gives us confidence in the results, even 
though we are aware of the fact that some assumptions are 
poorly constrained. The large majority of the stars have mass 
loss rates in a narrow range centered at ~3x10 _7 Mq yr" 1 , 
and it appears that very few AGB carbon stars (<5%) lose mat- 
ter at a rate less than ^5x 10 _8 Mq yr" 1 . 

We find that the mass loss rate and the gas expansion veloc- 
ity are relatively well correlated, but the scatter is large enough 
that we may conclude that the mass loss mechanism is able to 
produce a wide range of mass loss rates for a given expansion 
velocity. The mass loss rate is also well correlated with the pul- 
sational period of the star, correlated with the stellar luminos- 
ity, and there is a weak trend with the stellar effective tempera- 
ture, in the sense that the cooler stars tend to have higher mass 
loss rates. Also the gas expansion velocity is positively cor- 
related with the period and the luminosity. We conclude that 
the mass loss rate increases with increased regular pulsation 
and/or luminosity, and that the expansion velocity increases 
with mass loss rate (for low mass loss rates) and luminosity. 
The observed trends are all supporting the common concensus 
that these winds are driven by radiation pressure on dust grains, 
and that pulsation may play an important role. Somewhat sur- 
prising there appears to be no dependence on the stellar C/O- 
ratio. 

Our standard CSE model, assuming a single smooth ex- 
panding wind produced by a continuous mass loss, fails to re- 
produce the observational data for about 10% of the sample 
stars. Most notable among these are five stars with detached 
CSEs, presumably formed during a period of highly increased 
mass loss, and low present mass loss rates. We have found in- 
dications that the present mass loss rate is significantly higher 
for the star with the youngest dCSE, and that the shell mass 
may increase with shell age. Since our sample is reasonably 
complete, we can estimate that the time scale between mass 
ejections is about 10 5 years, if it is a repeatable phenomenon. 
An association with He-shell flashes is favoured. These objects 
have been the targets of a number of extensive, observational 
and theoretical studies. 

For some of our sample stars there exist enough observa- 
tional constraints to determine a combined dust parameter from 
the kinetic temperature structure. This result can be interpreted 
as a dust-to-gas mass ratio which increases with mass loss rate, 
but changes in the dust properties may also play a role. We also 
find that this means that the gas kinetic temperature in a carbon 
star CSE depends only weakly on the mass loss rate. 

The size of the CO envelope is an important parameter 
in the mass loss rate determination, at least for the low mass 
loss rate objects. We have used published radial CO( J=l— >0) 
brightness distributions, in combination with the radiative 
transfer code, to show that the CO photodissociation calcula- 



tion by Mamon et al. (1988) gives reasonably accurate results. 



In a few cases the observed radial brightness distributions are 
clearly different than the model results, suggesting deviations 
from a simple r" 2 density law, and hence time-variable mass 
loss. 
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We estimate that carbon stars, of the type studied here, re- 
turn on the order of 0.05 Mq yr _1 of gas to the ISM making 
them marginally important for the gas cycle in galaxies. More 
extreme carbon stars may contribute an order of magnitude 
more. However, they are probably not important as regards the 
origin of carbon. 
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Fig. 6. Multi-transition CO millimetre-wave line emission observed towards the bright carbon star W Ori. The observed spectra 
(histograms) have been overlayed with the model prediction (full line) using a mass loss rate of7xlO -8 M Q yr~ 1 . The transition, 
telescope used, and the corresponding beamsize, are indicated for each of the observations. The feature at ^5kms _1 may be 
of interstellar origin. The double-peaked profile in the CO( J=l— >0) model spectra is due to maser action and not due to spatial 
resolution effects. Cooling panel: The full line represents adiabatic cooling; the dotted line gives the H2 cooling; and the dashed 
line is CO cooling (see text for details). Temperature/optical-depth panel: The dotted line shows the kinetic gas temperature 
derived from the energy balance equation. The full line gives the excitation temperature of the CO( J=2— >1) transition, and a 
dash-dot line indicates a negative excitation temperature, i.e., maser action. The dashed lines give the tangential optical depth, 
T( an , of the CO( J=2— >1) transition. 
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Fig. 7. x 2 maps showing the sensitivity of the model to the adjustable parameters, i.e., the mass loss rate (A/) and the /i-parameter. 
Contour levels are drawn at Xmin+(l-0> 2.3, 4.6) with the middle contour (thick line) indicating the 68% confidence level. The 
best-fit model, as reported in Tab. ^, is indicated by a star and its corresponding % 2 and the number of observational constraints 
are also shown. 
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Fig. 9. Observed averaged radial brightness distributions overlayed by the results from the model, using the envelope sizes 
obtained from the chemical model by Mamon et al. (1988), for the parameters presented in Tab. [}| The circular beam used in the 
radiative transfer calculation is indicated by the dot-dashed line. Observational errors include both statistical errors and a 20% 
calibration uncertainty. 
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Fig. 10. Histograms showing the distribution of the mass loss rates for the whole sample (upper left panel), as well as subdivided 
into the different variability types (excluding LP And). Stars with known detached shells have not been included in these plots. 
In the upper left panel the "cross-hatched" area indicates the mass loss rate distribution for the stars within 500 pc of the Sun, 
i.e., the complete sample. 
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Fig. 12. The derived mass loss rate plotted against the luminosity (L), pulsational period (P), effective temperature (T c g), and 
the photospheric C/O-ratio of the star. The measured expansion velocity is plotted gainst the period and the luminosity. Sources 
with known detached shells are indicated with triangles (open triangles represent dCSEs, while filled ones represent aCSEs). 
Also shown are the power law index £, assuming the ordinate to scale with the abscissa, as well as the correlation coefficient r 
(the detached shell sources are not included in the fit). 



